Optimal Collision Energies and Bioinformatics Tools for Efficient Bottom-up Sequence Validation of Monoclonal Antibodies.

Rigorous validation of amino acid sequence is fundamental in the characterization of original and biosimilar protein biopharmaceuticals. Widely accepted workflows are based on bottom-up mass spectrometry, and they often require multiple techniques and significant manual work. Here, we demonstrate that optimization of a set of tandem mass spectroscopy (MS/MS) collision energies and automated combination of all available information in the measurements can increase the sequence validated by one technique close to the inherent limits. We created a software (called "Serac") that consumes results of the Mascot database search engine and identifies the amino acids validated by bottom-up MS/MS experiments using the most rigorous, industrially acceptable definition of sequence coverage (we term this "confirmed sequence coverage"). The software can combine spectra at the level of amino acids or fragment ions to exploit complementarity, provides full transparency to justify validation, and reduces manual effort. With its help, we investigated collision energy dependence of confirmed sequence coverage of individual peptides and full proteins on trypsin-digested monoclonal antibody samples (rituximab and trastuzumab). We found the energy dependence to be modest, but we demonstrated the benefit of using spectra taken at multiple energies. We describe a workflow based on 2-3 LC-MS/MS runs, carefully selected collision energies, and a fragment ion level combination, which yields ∼85% confirmed sequence coverage, 25%-30% above that from a basic proteomics protocol. Further increase can mainly be expected from alternative digestion enzymes or fragmentation techniques, which can be seamlessly integrated to the processing, thereby allowing effortless validation of full sequences.

Structural and computational basis for potent inhibition of glutamate carboxypeptidase II by carbamate-based inhibitors.

A series of carbamate-based inhibitors of glutamate carboxypeptidase II (GCPII) were designed and synthesized using ZJ-43, N-[[[(1S)-1-carboxy-3-methylbutyl]amino]carbonyl]-l-glutamic acid, as a molecular template in order to better understand the impact of replacing one of the two nitrogen atoms in the urea-based GCPII inhibitor with an oxygen atom. Compound 7 containing a C-terminal 2-oxypentanedioic acid was more potent than compound 5 containing a C-terminal glutamic acid (2-aminopentanedioic acid) despite GCPII's preference for peptides containing an N-terminal glutamate as substrates. Subsequent crystallographic analysis revealed that ZJ-43 and its two carbamate analogs 5 and 7 with the same (S,S)-stereochemical configuration adopt a nearly identical binding mode while (R,S)-carbamate analog 8 containing a d-leucine forms a less extensive hydrogen bonding network. QM and QM/MM calculations have identified no specific interactions in the GCPII active site that would distinguish ZJ-43 from compounds 5 and 7 and attributed the higher potency of ZJ-43 and compound 7 to the free energy changes associated with the transfer of the ligand from bulk solvent to the protein active site as a result of the lower ligand strain energy and solvation/desolvation energy. Our findings underscore a broader range of factors that need to be taken into account in predicting ligand-protein binding affinity. These insights should be of particular importance in future efforts to design and develop GCPII inhibitors for optimal inhibitory potency.

Selection of Collision Energies in Proteomics Mass Spectrometry Experiments for Best Peptide Identification: Study of Mascot Score Energy Dependence Reveals Double Optimum.

Collision energy is a key parameter determining the information content of beam-type collision induced dissociation tandem mass spectrometry (MS/MS) spectra, and its optimal choice largely affects successful peptide and protein identification in MS-based proteomics. For an MS/MS spectrum, quality of peptide match based on sequence database search, often characterized in terms of a single score, is a complex function of spectrum characteristics, and its collision energy dependence has remained largely unexplored. We carried out electrospray ionization-quadrupole-time of flight (ESI-Q-TOF)-MS/MS measurements on 2807 peptides from tryptic digests of HeLa and E. coli at 21 different collision energies. Agglomerative clustering of the resulting Mascot score versus energy curves revealed that only few of them display a single, well-defined maximum; rather, they feature either a broad plateau or two clear peaks. Nonlinear least-squares fitting of one or two Gaussian functions allowed the characteristic energies to be determined. We found that the double peaks and the plateaus in Mascot score can be associated with the different energy dependence of b- and y-type fragment ion intensities. We determined that the energies for optimum Mascot scores follow separate linear trends for the unimodal and bimodal cases with rather large residual variance even after differences in proton mobility are taken into account. This leaves room for experiment optimization and points to the possible influence of further factors beyond m/ z.

Pathways for Arene Oxidation in Non-Heme Diiron Enzymes: Lessons from Computational Studies on Benzoyl Coenzyme A Epoxidase.

Oxygenation of aromatic rings using O2 is catalyzed by several non-heme carboxylate-bridged diiron enzymes. In order to provide a general mechanistic description for these reactions, computational studies were carried out at the ONIOM(B3LYP/BP86/Amber) level on the non-heme diiron enzyme benzoyl coenzyme A epoxidase, BoxB. The calculations revealed four possible pathways for attacking the aromatic ring: (a) electrophilic (2e-) attack by a bis(µ-oxo)-diiron(IV) species (Q pathway); (b) electrophilic (2e-) attack via the σ* orbital of a µ-η2:η2-peroxo-diiron(III) intermediate (Pσ* pathway); (c) radical (1e-) attack via the π*-orbital of a superoxo-diiron(II,III) species (Pπ* pathway); (d) radical (1e-) attack of a partially quenched bis(µ-oxo)-diiron(IV) intermediate (Q' pathway). The results allowed earlier work of de Visser on olefin epoxidation by diiron complexes and QM-cluster studies of Liao and Siegbahn on BoxB to be put into a broader perspective. Parallels with epoxidation using organic peracids were also examined. Specifically for the BoxB enzyme, the Q pathway was found to be the most preferred, but the corresponding bis(µ-oxo)-diiron(IV) species is significantly destabilized and not expected to be directly observable. Epoxidation via the Pσ* pathway represents an energetically somewhat higher lying alternative; possible strategies for experimental discrimination are discussed. The selectivity toward epoxidation is shown to stem from a combination of inherent electronic properties of the thioacyl substituent and enzymatic constraints. Possible implications of the results for toluene monooxygenases are considered as well.

Mono- and binuclear non-heme iron chemistry from a theoretical perspective.

In this minireview, we provide an account of the current state-of-the-art developments in the area of mono- and binuclear non-heme enzymes (NHFe and NHFe2) and the smaller NHFe(2) synthetic models, mostly from a theoretical and computational perspective. The sheer complexity, and at the same time the beauty, of the NHFe(2) world represents a challenge for experimental as well as theoretical methods. We emphasize that the concerted progress on both theoretical and experimental side is a conditio sine qua non for future understanding, exploration and utilization of the NHFe(2) systems. After briefly discussing the current challenges and advances in the computational methodology, we review the recent spectroscopic and computational studies of NHFe(2) enzymatic and inorganic systems and highlight the correlations between various experimental data (spectroscopic, kinetic, thermodynamic, electrochemical) and computations. Throughout, we attempt to keep in mind the most fascinating and attractive phenomenon in the NHFe(2) chemistry, which is the fact that despite the strong oxidative power of many reactive intermediates, the NHFe(2) enzymes perform catalysis with high selectivity. We conclude with our personal viewpoint and hope that further developments in quantum chemistry and especially in the field of multireference wave function methods are needed to have a solid theoretical basis for the NHFe(2) studies, mostly by providing benchmarking and calibration of the computationally efficient and easy-to-use DFT methods.

Comparison of human glutamate carboxypeptidases II and III reveals their divergent substrate specificities.

UNLABELLED: Glutamate carboxypeptidase III (GCPIII) is best known as a homologue of glutamate carboxypeptidase II [GCPII; also known as prostate-specific membrane antigen (PSMA)], a protease involved in neurological disorders and overexpressed in a number of solid cancers. However, mouse GCPIII was recently shown to cleave ß-citrylglutamate (BCG), suggesting that these two closely related enzymes have distinct functions. To develop a tool to dissect, evaluate and quantify the activities of human GCPII and GCPIII, we analysed the catalytic efficiencies of these enzymes towards three physiological substrates. We observed a high efficiency of BCG cleavage by GCPIII but not GCPII. We also identified a strong modulation of GCPIII enzymatic activity by divalent cations, while we did not observe this effect for GCPII. Additionally, we used X-ray crystallography and computational modelling (quantum and molecular mechanical calculations) to describe the mechanism of BCG binding to the active sites of GCPII and GCPIII, respectively. Finally, we took advantage of the substantial differences in the enzymatic efficiencies of GCPII and GCPIII towards their substrates, using enzymatic assays for specific detection of these proteins in human tissues. Our findings suggest that GCPIII may not act merely as a complementary enzyme to GCPII, and it more likely possesses a specific physiological function related to BCG metabolism in the human body. DATABASE: The X-ray structure of GCPII Glu424Ala in complex with BCG has been deposited in the RCSB Protein Data Bank under accession code 5F09.

Transformation of Zwitterionic Pyridine Derivatives to a Spiro-Fused Ring System: Azoniabenzo[de]fluorine. Synthesis and Mechanistic Rationalization.

Reaction of aryl- and benzylsulfanopyridinium amidates bearing a methyl group in position 6 with 2 equiv of diphenylketene afforded a spiro-fused ring system: azoniabenzo[de]fluorine. By use of an excess amount of ketene, a distinct reaction was observed via which a 1H-pyrrolo[3,2-b]pyridin-2(3H)-one derivative was furnished. The structure of the tetracyclic spiro-fused ring system was unambiguously confirmed by X-ray diffraction, and its formation was rationalized by DFT calculations.

Reactivity of the binuclear non-heme iron active site of Δ9 desaturase studied by large-scale multireference ab initio calculations.

The results of density matrix renormalization group complete active space self-consistent field (DMRG-CASSCF) and second-order perturbation theory (DMRG-CASPT2) calculations are presented on various structural alternatives for the O-O and first C-H activating step of the catalytic cycle of the binuclear nonheme iron enzyme Δ(9) desaturase. This enzyme is capable of inserting a double bond into an alkyl chain by double hydrogen (H) atom abstraction using molecular O2. The reaction step studied here is presumably associated with the highest activation barrier along the full pathway; therefore, its quantitative assessment is of key importance to the understanding of the catalysis. The DMRG approach allows unprecedentedly large active spaces for the explicit correlation of electrons in the large part of the chemically important valence space, which is apparently conditio sine qua non for obtaining well-converged reaction energetics. The derived reaction mechanism involves protonation of the previously characterized 1,2-µ peroxy Fe(III)Fe(III) (P) intermediate to a 1,1-µ hydroperoxy species, which abstracts an H atom from the C10 site of the substrate. An Fe(IV)-oxo unit is generated concomitantly, supposedly capable of the second H atom abstraction from C9. In addition, several popular DFT functionals were compared to the computed DMRG-CASPT2 data. Notably, many of these show a preference for heterolytic C-H cleavage, erroneously predicting substrate hydroxylation. This study shows that, despite its limitations, DMRG-CASPT2 is a significant methodological advancement toward the accurate computational treatment of complex bioinorganic systems, such as those with the highly open-shell diiron active sites.

Reactivity models of hydrogen activation by frustrated Lewis pairs: synergistic electron transfers or polarization by electric field?

Two alternative qualitative reactivity models have recently been proposed to interpret the facile heterolytic cleavage of H2 by frustrated Lewis pairs (FLPs). Both models assume that the reaction takes place via reactive intermediates with preorganized acid/base partners; however, they differ in the mode of action of the active centers. In the electron transfer (ET) model, the hydrogen activation is associated with synergistic electron donation processes with the simultaneous involvement of active centers and the bridging hydrogen, showing similarity to transition-metal-based and other H2-activating systems. In contrast, the electric field (EF) model suggests that the heterolytic bond cleavage occurs as a result of polarization by the strong EF present in the cavity of the reactive intermediates. To assess the applicability of the two conceptually different mechanistic views, we examined the structural and electronic rearrangements as well as the EFs along the H2 splitting pathways for a representative set of reactions. The analysis reveals that electron donations developing already in the initial phase are general characteristics of all studied reactions, and the related ET model provides qualitative interpretation for the main features of the reaction pathways. On the other hand, several arguments have emerged that cast doubt on the relevance of EF effects as a conceptual basis in FLP-mediated hydrogen activation.

Hydrogen activation by frustrated lewis pairs: insights from computational studies.

Sterically encumbered Lewis acid-base pairs, the so-called frustrated Lewis pairs, can split dihydrogen heterolytically and act as transition metal free catalysts for the hydrogenation of unsaturated compounds. Here we review the results from our quantum chemical calculations aimed at the understanding of this remarkable class of reactions and we put them into the context of related works from other research groups. The thermodynamics of the H2 splitting reaction is discussed first; the role of acid-base properties, intramolecular cooperativity, and other factors is assessed, employing an energy partitioning scheme and also in the light of the latest experimental findings. The mechanism of hydrogen cleavage is then examined, and an overview about the applicability of our reactivity model involving synergistic electron transfers between H2 and preorganized Lewis acid/base centers is given. Finally, insights about catalytic cycles in FLP-mediated hydrogenations are summarized, pinpointing the diversity of the involved elementary steps and their possible sequences.

Coordination and bond activation in complexes of regioisomeric phenylpyridines with the nickel(II) chloride cation in the gas phase.

Electrospray ionization of dilute solutions of phenylpyridines (phpy) in the presence of nickel(II) chloride leads to gaseous ions of the type [Ni(phpy)(m)](2+) with m = 3-5 and [NiCl(phpy)(n)](+) with n = 1-3, which are characterized by various gas-phase experiments in combination with calculations using density functional theory. Of the regioisomeric phpy's, 2-phpy behaves drastically different compared to 3- and 4-phpy. Ion mobility mass spectrometry allows a differentiation of the gaseous ions and an elucidation of characteristic properties of the metal complexes. For 2-phpy, C-H bond activation in the [NiCl(phpy)(2)](+) complex is significant, whereas this route is almost suppressed for the corresponding complexes of 3- and 4-phpy and only occurs at elevated energies.

Identification and interconversion of diastereomeric oligo-Tröger bases probed by ion mobility mass spectrometry.

Oligo-Tröger bases are auspicious scaffolds of molecular engineering, which motivates studies on the mechanism of their interconversion and on the facile determination of the relative configuration of their diastereoisomers. Protonated, sodiated, and argentated species of those compounds were therefore studied via ion-mobility mass spectrometry (IM-MS), allowing differentiation on the basis of the shapes of the ions. First, the isomerization was confirmed to be acid-catalyzed as it takes place readily in the case of protonated Tröger bases, whereas the metallated bases are configurationally stable. Second, the corrected arrival times of the various isomers of the cationized bases were found to show distinct differences in IM-MS, and their excellent correlation with the cross sections obtained from quantum chemical calculations paves the way toward the easy identification of diastereoisomers.

Theoretical calculations of physico-chemical and spectroscopic properties of bioinorganic systems: current limits and perspectives.

In the last decade, we have witnessed substantial progress in the development of quantum chemical methodologies. Simultaneously, robust solvation models and various combined quantum and molecular mechanical (QM/MM) approaches have become an integral part of quantum chemical programs. Along with the steady growth of computer power and, more importantly, the dramatic increase of the computer performance to price ratio, this has led to a situation where computational chemistry, when exercised with the proper amount of diligence and expertise, reproduces, predicts, and complements the experimental data. In this perspective, we review some of the latest achievements in the field of theoretical (quantum) bioinorganic chemistry, concentrating mostly on accurate calculations of the spectroscopic and physico-chemical properties of open-shell bioinorganic systems by wave-function (ab initio) and DFT methods. In our opinion, the one-to-one mapping between the calculated properties and individual molecular structures represents a major advantage of quantum chemical modelling since this type of information is very difficult to obtain experimentally. Once (and only once) the physico-chemical, thermodynamic and spectroscopic properties of complex bioinorganic systems are quantitatively reproduced by theoretical calculations may we consider the outcome of theoretical modelling, such as reaction profiles and the various decompositions of the calculated parameters into individual spatial or physical contributions, to be reliable. In an ideal situation, agreement between theory and experiment may imply that the practical problem at hand, such as the reaction mechanism of the studied metalloprotein, can be considered as essentially solved.

Structural and spectroscopic properties of the peroxodiferric intermediate of Ricinus communis soluble Δ9 desaturase.

Large-scale quantum and molecular mechanical methods (QM/MM) and QM calculations were carried out on the soluble Δ(9) desaturase (Δ(9)D) to investigate various structural models of the spectroscopically defined peroxodiferric (P) intermediate. This allowed us to formulate a consistent mechanistic picture for the initial stages of the reaction mechanism of Δ(9)D, an important diferrous nonheme iron enzyme that cleaves the C-H bonds in alkane chains resulting in the highly specific insertion of double bonds. The methods (density functional theory (DFT), time-dependent DFT (TD-DFT), QM(DFT)/MM, and TD-DFT with electrostatic embedding) were benchmarked by demonstrating that the known spectroscopic effects and structural perturbation caused by substrate binding to diferrous Δ(9)D can be qualitatively reproduced. We show that structural models whose spectroscopic (absorption, circular dichroism (CD), vibrational and Mössbauer) characteristics correlate best with experimental data for the P intermediate correspond to the µ-1,2-O(2)(2-) binding mode. Coordination of Glu196 to one of the iron centers (Fe(B)) is demonstrated to be flexible, with the monodentate binding providing better agreement with spectroscopic data, and the bidentate structure being slightly favored energetically (1-10 kJ mol(-1)). Further possible structures, containing an additional proton or water molecule are also evaluated in connection with the possible activation of the P intermediate. Specifically, we suggest that protonation of the peroxide moiety, possibly preceded by water binding in the Fe(A) coordination sphere, could be responsible for the conversion of the P intermediate in Δ(9)D into a form capable of hydrogen abstraction. Finally, results are compared with recent findings on the related ribonucleotide reductase and toluene/methane monooxygenase enzymes.

Curvature correction for microiterative optimizations with QM/MM electronic embedding.

One of the most common methods to treat the electrostatic effect of the environment in QM/MM calculations is to include the MM atoms as point charges in the QM Hamiltonian. In this case, a microiterative geometry optimization ignoring the QM contributions to the forces in the relaxation of the environment cannot yield exact stationary points. One solution that has been suggested in the literature is based on using a constant additive correction to the MM gradient during the microiterations, determined in the preceding macroiteration. Here, we analyze the convergence properties of the gradient correction method and point out that a smooth relaxation is not ensured if the curvature of the approximate, MM-based description of the potential energy surface of the environment is too small in comparison with the exact one. We suggest a computationally cheap second-order correction that uses an estimated Hessian from the Davidon-Fletcher-Powell method to tackle the problems caused by the too small curvature. Test calculations on four metalloenzymatic systems (∼100 QM atoms, ∼2000 relaxed MM atoms, ∼20,000 atoms in total) show that our approach efficiently restores the convergence where gradient correction alone would lead to oscillations.

Rearrangement of aryl- and benzylthiopyridinium imides with participation of a methyl substituent.

6-Methyl substituted 2-aryl- and 2-benzylthiopyridinium N-imides reacted with an excess of isocyanates to give N,N-disubstituted exocyclic1H-imidazo[4,5-b]pyridin-2(3H)-ones. The products easily underwent spontaneous [1,5] hydrogen shift to provide the heteroaromatic imidazopyridinone isomers. The transformation implied the initial formation of [1,2,4]triazolo[2,3-a]pyridinium salt, followed by deprotonation and carbamoylation of the methylene moiety, and, finally, a rearrangement following a [1,3] sigmatropic pattern. Mechanistic considerations suggest and some experimental findings reveal the nonconcerted two-step mechanism of the ring transformation step.

On the mechanism of the copper-mediated C-S bond formation in the intramolecular disproportionation of imine disulfides.

The mechanism of the copper-mediated disproportionation of aromatic imine disulfides to benzothiazoles in the gas phase is investigated by experimental and theoretical methods. Application of infrared multiphoton dissociation and hydrogen/deuterium exchange experiments combined with density functional theory (DFT) calculations of the relevant molecular structures and the associated infrared spectra allows the identification of the observed ionic intermediates. The theoretical investigation of the possible reaction pathways supported by collision-induced dissociation experiments provides a consistent mechanistic picture of the reaction catalyzed by a single copper(I) ion. Activation of the substrate proceeds via homolytic sulfur-sulfur bond cleavage, yielding metal complexes in the formal +3 oxidation state; carbon-sulfur coupling and hydrogen-atom transfer complete the transformation to the products. Exploratory studies demonstrate that in the gas phase, the disproportionation of the imine disulfide can also be mediated by other metal ions via different either homo- or heterolytic mechanisms without involving high-valent intermediates.

The solid state structure and reactivity of NbCl(5) x (N,N'-dicyclohexylurea) in solution: evidence for co-ordinated urea dehydration to the relevant carbodiimide.

NbCl(5) x (N,N'-dicyclohexylurea) 1a owns a distorted octahedral structure due to intramolecular NH...Cl bonding. The unit cell contains four units which are intermolecularly NH...Cl and NH...N bonded. An extended intramolecular network of H-bonding (N-H...Cl, CH...Cl, CH...N) causes the 3D self assembling of the units. Upon addition of base, the HCl release from 1a is observed with the transfer to Nb of the O-atom of the carbonylic function of the starting urea which is converted into the relevant carbodiimide CyN=C=NCy 4. The latter is quantitatively released by adding an excess of NEt(3) at 308 K (py and DBU are less efficient) with formation of the known NbOCl(3)(NEt(3))(2), isolated in quantitative yield. Increasing the temperature leads to a loss in selectivity as the formed DCC undergoes further reactions. At 350 K, the isocyanate CyN=C=O has been isolated in 60% yield besides a mixture of Nb-complexes. DFT calculations have been coupled to IR and NMR experiments for characterizing possible reaction intermediates and the behaviour of 1a. Several other MCl(x) species (ScCl(3), YCl(3), LaCl(3), TiCl(4), TaCl(5), AlCl(3), SnCl(4)) have been shown to be able to co-ordinate DCU but not all of them promote the conversion of urea into DCC.